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Abstract 

It is shown that instability of homogeneous precession is caused by 
combined effect of anisotropy of spin wave velocities and dipole inter- 
action. In the principal order on the ratio of the Leggett frequency 
to the Larmor frequency the increments of growth of spin wave am- 
plitudes are found. The magnitude of the maximum increment for all 
deviation angles of spin from its equilibrium orientation is calculated. 
The estimation is made of the minimum temperature down to which 
the precession is stable. 

1. In a contrast to the superfluid A-phase, in the B-phase long wavelength 
perturbations do not destroy homogeneous precession of magnetization. Pre- 
cession is stable with respect to such perturbation if the initial deviation angle 
of magnetization f3 > 6q , where 8q = arccos(— 1/4) and marginally stable if 
(3 < 9 . However, at low temperatures fast decay of homogeneous precession 
is observed in both cases. Such decay was first observed in experiments [U, [2] 
and is referred as catastrophic relaxation. Earlier we considered a process 
of parametric excitation of spin waves with finite wave vectors by the pre- 
cession of magnetization as a possible origin of the decay j3]. In the theory 
of magnetics this effect is known as the Suhl instability pi]. The distinctive 
property of the instability in 3 He-B is that it occurs at the precession with 
large tipping angles (3 ~ 100° and also that excitation of different types of 
spin waves is possible. In Ref. [3] a scheme was proposed that takes into 
account both mentioned peculiarities and increments of instability for each 
of three types of spin waves were found. The dependence of the obtained 
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increments on magnetic field does not agree with that experimentally ob- 
served [5|. The magnitudes of the increments were also overvalued. The 
disagreement with the experiment was caused by our technical mistake. As 
a result of this mistake it turned out, that the instability of the precession 
can be induced by anisotropy of spin wave velocities in 3 He-B alone. The 
further analysis has shown that in order to provide coupling between spin 
waves and precession the dipole interaction has to be taken into account. It 
turned out also that considerable contribution to the instability comes from 
joint resonances when spin waves belonging to different branches of spectra 
are excited simultaneously. 

In the present paper a revised version of the theory of parametric in- 
stability of homogeneous precession of spin in 3 He-B is presented. At the 
same time the foregoing disagreement with the experiment is resolved. A 
comparison is also made with the results of Ref . (6j [7] , where contribution of 
the boundaries to the development of parametric instability is considered. 

2. Following the procedure of Ref. [3] let us parameterize orientation of 
the order parameter of 3 He-B, which is a rotation matrix , by the Euler 
angles a , ft , 7 ( z -axis is oriented opposite to the direction of the d.c. mag- 
netic field Ho). Actually, it is more convenient to use the sum $ = a + 7 
instead of the angle 7 . Canonically conjugated momenta to these coor- 
dinates are correspondingly the following combinations of spin projections 
P = S z — , Sp , , where S z — is projection of spin onto z-axis, - 
its projection onto ( = Rz and Sp — is projection on the line of nodes (see 
for example [8]). Equations of motion are hamiltonian with respect to the 
stated pairs of variables with the Hamiltonian 

(1) 

Here ujl - - is the Larmor frequency, corresponding to the d.c. magnetic 
field, F v — the gradient energy, Ud(a,ft, <3>) — the dipole energy, which is 
of the order of the squared Leggett frequency in the 3 ife-B: VL 2 B . We choose 
the units of measurements so that the magnetic susceptibility of 3 He — B - 
X , and the gyromagnetic ratio for nuclei of 3 He — g are equal to unity. In 
this units spin has a dimensionality of frequency and energy — of a squared 
frequency correspondingly. In the standard setting of NMR experiments in 
the regions distant from the walls of the cell spin precesses in the so called 
Leggett configuration, when the "orbital vector" U = —R^s^ [6] is parallel 
to the magnetic field. In this case Ud(oc, ft, $) does not depend on the angle 
a and precession is described by the stationary solution of equations of spin 
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dynamics which do not contain oscillating terms: 

a = —uj p t + a , 7 = ujpt + $^°- ) — a , 
P (0) = u p {cos(3 - 1), Sf ] = 0, Si 0) = Up cos /3, (2) 

where u; p - the frequency of precession. If (3 > 0q then = , and if 

cos$ (0) = (^-cos/? (0) )/(i + cos/3 (0) ). (3) 

Explicit time dependence of the stationary solution (j2J) can be excluded if one 
transfers to the variable if> = a + u p t and uses new Hamiltonian H = H + 
6c>pP , then ^ = . To find the spectra of excitation against the background 
of precession we linearize the equations of motion on small deviations from 
the stationary solution (EJ E|): 5ip(r,t) — ip — vp^ , etc. For the sake of 
convenience the following combinations of the mentioned deviations are used: 

5P +(1- cos /3)5S< 

v = sin pdip, v = -, 

ul sin p 

£ = <5$-(l-cos/3)5V>, a = 5S c /u L , (4) 
( = 5Sp/u L , 7] = -5(3. 

Expression for the gradient energy of 3 He — B contains two coefficients that 
can be written as velocities of two types of spin waves cm and c± . In what 
follows units of length and time are chosen so that ul = 1 and c\\ = 1 . 
Without loss of generality one can assume that variables change only in y 
and z directions, then the time-independent part of the gradient energy has 
the form: 

Fvst = ~[(1 - AOW + Vl) + "I + V 2 Z + 4 + (1 - ^)e% (5) 

where \i — \ — c\/c 2 ^ - is the anisotropy of spin wave velocities. Parameter 
H will be considered as a small one, in fact fi fsl/4 [9]. Furthermore, Fy 
also contains term oscillating with the frequency of precession: 

Fvosd = -V[(v,y £ ,z + V,z£, y ) cos Upt + (v, y e iZ + v, z e, y ) sxnu p t) (6) 
and with the doubled frequency of precession 

Pvosc2 = ~ [{v 2 y - v 2 y ) cos 2u p t + 2v t y-q >y sin 2u p t] . (7) 

In zero order approximation on small parameters [i and (Qb/^l) 2 the equa- 
tions for deviations have hamiltonian form with the Hamiltonian: 

h = \[& + V) 2 + C 2 + ° 2 + (W) 2 + (Vr?) 2 + (Ve) 2 ] (8) 
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and with respect to the pairs of canonically conjugated variables (e, a) ; 
(i/, $) ; ((, if) . In each pair the first variable is coordinate and the second is 
momentum. Equations of motion for the pair (e, a) have the form: 

at = (9) 

It is convenient to rewrite them in a vectorial form: 

dXi 



dt 



Af x Xi. (10) 



Solutions of the system ([9]) have the form of plane waves ef exp[i(kr =)= uit)] 
with the dispersion law 

u x = k. (11) 

Here ef are right eigenvectors of matrix M\ corresponding to eigenvalues 
=RWi : 

ef=( 1 V (12) 

One needs left eigenvectors fj* 1 of the same matrix to be able to make pro- 
jections. They can be normalized so that the following conditions are met: 

(tf,ef) = l; (tf,eT) = 0. (13) 

Here the scalar product is defined as 

(a,b) = £ a* n b n . (14) 

n=l,2 

As a result 

^=2^(wi;=F<)- ( 15 ) 

In a similar way the following dispersion laws for two transverse modes are 
obtained: 

«* = J\ + * l -\, (16) 



<"3 = ^ + * 2 + ^ (17) 

At /c — > co> 2 ~ A^ 2 , i.e. it is a gapless mode arising from the degeneracy 
of precession with respect to a . Another mode has a gap u 3 = lol at 
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k — > , it passes into nutations. In the coordinates ip, $, rj, ( the following 
right eigenvectors correspond to the frequencies ±uj2 '■ 



and to the frequencies ±003 : 



( 1 \ 

±i 
V ^2 J 



I 1 \ 

=R^2 
V -W3 / 



(19) 



Left eigenvectors are correspondingly 

1 



L 2 



2u 
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and 



r 3 



2cj 



cj 3 ; =Fi; ±iu 2 ; 1) , 



-cj 2 ; ±i; ±iuj 3 ] 1) 



(20) 



(21) 
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where u 2 3 = u> 2 + uj 3 . 

3. Time-dependent corrections to the Hamiltonian (jSj) can provide cre- 
ation and mutual transformation of excitations. At l||Ho the dipole energy 
does not contain time-dependent terms and the gradient energy in the first 
approximation on \x contains oscillating terms (EI), (CO). Taking the oscil- 
lating terms into account one can write equations of motion for deviations 
e, a, v, 1?, (, 7] combined in a six-component vector-column in a form: 



^ = (Mo + V(tj) X, 



(22) 



where all time-dependent terms are collected in V(t) . The sum of ([6]) and 



yields V(t) = £ 



n=l,2 



Following pro- 



W n exp(—inu p t) + W* exp(inujpt) 
cedure of time-dependent perturbation theory let us seek for a solution of the 
system f[2"2l in a form of expansion in eigenvectors of matrix M 



X(r, t) = y^{ a+ k-(^) e+ k- ex P(^kr — iujjt) + a.j^(t)e.j^exp(«kr + iujjt)} (23) 



supplemented with 



Eigenvectors are the mentioned above e\ , c 2 , c 3 
zeroes up to six component. Substitution of (1231 into ( 1221) and separation of 
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equations by k yields: 

I^^i^H^) + «7k e 7k ex P(^0} = 

3 

^(W,, exp(— intjjpt) + IV"^ exp(mto' p t)){a + j 4 .e + j 4 . exp(— iu/jt) + a _ j [ e _ ] t exp(ia; :) t)}(24) 

Multiplying two sides of equation (l2~4"j) by exp(iu;/t) and omitting common 
index k one has: 

aj = ^{ajffi, W n ej) exp[i(ui + ujj - nu p )t] + 

af(ff, W*e+) exp[i(u;, - + nu p )t}}. (25) 

Application of the procedure assumes that coefficients are weakly 

varying at time ~ l/u) e - Averaging of equation ( 1251) shows that non- 
trivial correlations between different a * arise only nearby the resonances 

C0j(k) — ui(k) = nujp and u>i(k) + k) = ncj p . The second resonance cor- 
responds to creation of quasiparticles from "vacuum". In accordance with 
equality ^(k) = 0Jj{— k) the sign of one of the momenta is changed. In 
what follows it is assumed that the resonance condition is fulfilled exactly 
then the obtained increment of instability is a maximum one. If at given k 
the resonance condition is fulfilled only for two states I and j , then 

a+ = {ft,Wej)aj. (26) 

The same argument for aj yields: 

aJ = (f7,W*e+)a+. (27) 

The system of equations ([211 . (1271) has solutions ~ exp(±At) , where A is 
determined by A 2 = (fj~, W / e^)(fJ f , Wej)* . Thus the problem of determina- 
tion of the increment of instability reduces to calculation of the elements 
of matrix (fj,Wef) between the states satisfying the resonance condi- 
tions. Particularly, if I = j , then 2dj/(k) = nu p and A 2 = |(ft, lye;")! 2 . 
This is the simplest case of parametric resonance [10J. Explicit expres- 
sions for W± and W 2 are found using equalities (jSJ) and (J7j). There 
are four non-zero elements of matrix W\\ (Wi) av = (Wi)& e = ifik y k z and 
{Wi) av = (Wi)^ £ = fik y k z , and there are another four finite elements of matrix 
W 2 : (W 2 )^=-(W 2 )< :v = fxkl/2 and (W 2 )^ = {W 2 ) Cv = . Here k y , k z 

are the components of the wave vector satisfying the resonance conditions. 
Direct verification shows that the matrix element (fj~, Wief) is not equal to 
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zero only if j — 1,1 — 3. The resonance condition u>i(k) + k) = ujl 
is fulfilled only at k = , but at this value of k the matrix element turns 
into zero due to the factor k y k z . For the matrix (f^ ,W2&i) non-zero el- 
ement corresponds to j = I = 3 . In this case the resonance condition 
cj 3 (k) + cj 3 (— k) = 2u>l is also fulfilled only at k = and the corresponding 
matrix element is equal to zero due to the factor ky . Thus if the dipole en- 
ergy is neglected then the anisotropy of spin wave velocities does not provide 
the coupling between precession and spin waves. 

4. In order to take the dipole energy into account it is need to add to 
the RHS of equation (1221) the matrix of dipole torque Nd and to find new 
energies of excitations and the corresponding eigenvectors. The added terms 
in the equation ([2"2l are small in comparison with the elements of the matrix 
Mq as (Qb/^l) 2 i which is in typical experimental conditions of the order of 
~ lCT 1 — 10~ 2 . In the first order approximation on the mentioned parameter 
the elements of matrix (ft, W n ef) are given by the expressions: 

(ft, W- n ef) = i ]T ($,N D e m )(f m ,W n er) + 

i E " ($,W n e m )(t m ,N D ei). (28) 



U0[_ — uo. 



in 



Non-zero elements are obtained at the following resonances: 

wi(Jfe) +Wi(-Jfe) =u p , k = 1/2, (3< 6 , (29) 

ui(k) +uj 2 {-k) = u p , k = 2/3, (30) 

uj^k) +u 3 (-k) = 2u p , k = 2/3, (31) 

u 2 {k) + u 3 (-k) = 2u p , k = V3/2. (32) 

The corresponding increments are of the order of [iQ 2 b /ljl . At A; =1/2 
( , v /(l-cos(/3))(l + 4cos(/3)) ^ 

1 /3 > 6 , 

where 5 is the angle between the direction of the wave vector k and the 
direction of the magnetic field. At k=2/3 there are two resonances and one 
needs to consider the system of equations for thee amplitudes in order to find 
the increment. In this case the increment is given by the expression 



A(12,13) - \Afl2) + A(i 3 ), (34) 
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Figure 1: 

Dependence of the coefficient a in formule flHHj) on tipping angle (3 



where 



A 



(12) 



. /nrx VlO(l-cos(/3))Q| a „ 

//sin(25)- — ^— — ^IZZ—B, /3<0 O , 



/i sin (25) 



25 cj l 
|4cos 2 (/3) + 31cos(/3) + 15|fi| 



H13) 



[ism 2 (5) 



15V10 w L 
10^(1 -cos(/3))(l + 4cos(/?)) ft 



/3>0 O , 



25 



• 2 , n v / 10si n(/j)(l + 4cos(/?))Oj ? 
-/ism (o) 



225 



0<0o, 
P>0o, 



The increment for the resonance at k = V3/2 is 



A 



(23) 



H sin 2 (5) 



3(1 - cos(/3)) ft| 



20 



. 2 ^J4cos 2 (/?) + 31cos(/3) + 15|fi| 
/ism 2 (£) J - 



P<0o, 

, (3>e . 



(35) 



(36) 



(37) 



(3? 



40 u L 
The maximum increment for each value of (3 can be written as 

>^max((3) = a((3)fl — , 

where the dependence of the coefficient a on tipping angle is shown at Fig.l. 

At finite temperatures the damping of spin waves has to be taken into 
account. The instability sets up if the increment of growth of waves, which 
satisfy the resonance condition, exceeds the decrement of damping. As before 
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[3] for estimation of the temperature of catastrophic relaxation T cat it will 
be assumed here that the principal mechanism of dissipation is spin diffusion. 
The minimum temperature, down to which the precession is stable, is found 
from the equation: 

£^-X (39) 

where D{T) is the coefficient of spin diffusion. At temperatures in question 
T < 0.4T C , the increment weakly depends on temperature and its value can 
be taken at T = . The LHS of ( 1391) strongly depends on temperature owing 



to spin diffusion which behaves as D(T) ~ J T/A exp(— A/T) at T — > . 

The obtained increments originate from the coupling between precession 
and spin waves in bulk helium. In Ref. E] the increment arising from 
enhancement of the coupling in regions adjacent to the walls was found. 
Because of boundary conditions precession on the walls goes on in the con- 
figuration different from the Leggett one. In this case oscillating terms in the 
dipole energy appear without taking the anisotropy of spin waves velocities 
into account. The local coupling in regions adjacent to the walls is on the 
order of greater than that obtained for the bulk helium. In Ref. [HI E] the 
result of calculation of increment of growth and T cat for the angle (3 = 90° is 
represented. The bulk contribution to the increment for the same conditions 
estimated with the use of the above formulae gives approximately the same 
value as the surface contribution. It is impossible to separate surface and 
bulk contributions to the increment by their dependence on magnetic field 
because both contributions are proportional to VL 2 b /ujl . However, it has to 
be mentioned, that the surface contribution depends on ratio of volume ad- 
jacent to the surfaces to the total volume of helium. For angles (3 < 104° 
there is no characteristic length whereon penetrates the perturbation effect 
of the walls on the precession. Thus the regions adjacent to the walls oc- 
cupies considerable part of the volume. The most part of the data about 
catastrophic relaxation is obtained in experiments with the homogeneously 
precessing domain. In this case magnetization precesses with angles slightly 
above 8q rs 104° and the frequency of precession is shifted from the Larmor 
frequency. Then the effect of the walls is limited by the "coherence length" 
£ = cj_/ Jujp(u p — ujl) [HIE]- For the typical experimental conditions in 
the most part of the precessing domain £ ~ cm and the regions adjacent 
to the walls occupies only a small part of the total volume. The increment 
is determined by the bulk resonance o; 2 (A;) + u 3 (—k) = 2u p . The available 
experimental data for diffusion coefficient p3] pertains to the temperatures 
T > 0, 4T C . In the experiments at pressure P ~ 20 bar and magnetic field 
H ~ 142 Oe T cat ~ 0.42T C . After substitution of the diffusion coefficient 
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at this temperature D = 0.04 cm 2 /c and the values of other parameters 
Q = 2tt ■ 250 kHz,C|| = 1.6 ■ 10 3 cm/c, \i = 1/4 into the formula (JMD we 
obtain for the LHS of equation ( !39l) the value 5 • 10 4 1/c, and for the RHS 
- 4 -10 4 1/c. Taking into account that parameters \i and VL 2 b /uj 2 l are not 
too small the obtained agreement can be regarded as a satisfactory for this 
example. 

5. In conclusion the Suhl instability limits from below the interval of tem- 
peratures where coherent precession in 3 He — B can exist. The obtained here 
low temperature limit of stability of precession is caused by the interaction 
of precession with spin waves in the bulk helium. 

Lowering of the limit temperature can be achieved by using higher mag- 
netic fields as it was demonstrated in experiments Such tendency agrees 
with formula f[3"9l . On the one hand increasing of magnetic field decreases 
the increment of instability and on the other hand it increases spin waves 
damping. However, due to the exponential dependence of the coefficient of 
diffusion on temperature the effect of magnetic field on T cat becomes weaker 
when temperature decreases. Detailed comparison with the results of [H] will 
be done in a full-length publication. 

We thank to V.V. Dmitriev for useful discussion. This research was sup- 
ported by RFBR and Ministry of Education and Science of Russian Federa- 
tion. 
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